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ABSTRACT 

Empirical  Orthogonal  Function  (EOFI  rsprassntatioD  of  the 
synoptic  forcing  is  canbined  with  past  meridional  and  zonal 
dis placements  (persistence)  to  forecast  tropical  cyclone 
rscurvature  at  36,  54  and  72  h.  Recurvature  is  defined 
following  Leftwich  (1978,  1932):  Recirvature  is  a  net 
displacement  northward  of  SIS*^  luring  the  forecast  period  or 
the  attaining  of  northeastward  motion  daring  the  12  hours 
prior  -0  the  end  of  the  forecast  period.  The  combination  of 
persistence  and  EOF  coefficients  consistently  out-prrf ormed 
the  individual  methods  for  forecasting  recurvature  in  both 
the  dependent  and  independent  data  samples.  Evaluatioji 
scores  were  the  Brier  p-score,  the  Heiike  skill  score  and 
percent    correctly    forecast. 

Seven  forecast  aids  used  by  Joint  Typtioon  Warning  Cen^-r 
(JTWC)  Guam,  were  tested  f^r  ability  tD  forecast  -ropical 
cyclone  recurvature  at  2U,  48  and  72  h.  Using  the  skill 
scores  described  abcve  for  tie  1979-31  typhoon  seasons, 
performance  indicates  large  annual  varia~ions  of  skill 
between   -che   models. 
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I.        INTSDDOCTION 

Tropical  cyclone  recurvature  is  one  of  the  most 
difficult  events  to  forecast  for  even  tha  most  experienced 
forecaster.  The  ability  to  forecast  recurvature  correctly 
could  annually  save  both  the  military  and  civillian 
communities  millions  cf  dollars,  thousands  of  man  hours  and 
the    lives      of   many   victims.  There    are   few      reliable    rules 

that  the  tropical  cyclone  forecaster  can  follow  for 
forecasting  recurvature  cr  non-recurvature  in  either  the 
near  and  the  extended  time  frames.  Riehl  and  Shafer  (1944) 
were  among  the  first  to  summarize  guidelines  for  forecasting 
recuTvature  cf  Atlantic  tropical  cyclones.  Riehl  and  Shafer 
defined  recurvature  as  "any  major  alteration  of  a  storm 
track"  which  is  quite  non-descri pt.  RiahL  and  Shafer  listed 
several  factors  affecting  tropical  cyclone  motion,  including 
the  effects  of  the  upper-level  winds  and  the  interaction 
with  troughs  in  the  westerly  flow  to  -he  west  of  the 
tropical  cyclone.  They  also  describe!  various  synoptic 
situations  to  assist  the  forecaster  in  predicTing  ihe  storm 
path.  The  Joint  Typhoon  Warning  Center  (JPrtC)  Guam  defines 
recurvarure   as 
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"The  turning      of  a    tropical      cyclone   frDm    an      initial   path 
toward   the   west   and  northwest    to  the   north    then   northeast." 

which   is   also    quite   non-descr ipt . 

More  recent  studies  by  George  and  Gray  (1977),  Chan,  et 
al     (1980)       and   Shaffer    (1982)  have    shown    that    the   tropical 

cyclone  motion  is  largely  dependent  upDn  the  surrounding 
steering    flow    represented   by   height/temperature   fields. 

The  purpose  of  this  study  is  to  use  a  simple 
representation  of  the  synoptic  forcing  by  Empirical 
Orthogonal     Functions      (EOF's)  plus        persistence-related 

variables  to  predict  tropical  cyclone  recurvature  or 
non-recurvature.  The  time  increments  of  36,  54  and  72  h 
for  forecasting  recurvature  were  selected  for  several 
reasons.  First,  a  minimum  tiae  period  of  36  h  was  thought 
to  be  of  greatest  value  and  to  provide  the  most  confident 
forecast  for  the  shore  establishment,  air  wing  squadrons  and 
ships  presently  in  port.  Secondly,  The  apper-air  soundings 
are  transmitted  to  Fleet  Numerical  Oceanography  Center 
(FNOC)  by  approximately  3  h  after  00/12  GMT.  These 
soundings  are  used  tc  generate  the  12-hDurly  forecasts  of 
the  500  mb  heigh-s  (in  this  case).  These  generated  fields 
would   be      available    tc     Jcint    Typhoon      viarning    Center     (JTWC) 
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Guam  prior  to  the  next  warning  at  approximately  6  h 
after  00/12  GMT»  The  computer  program  to  decorapos-e  the 
D-value  fields  (deviations  from  the  standard  atomsphere 
gsopotential)  into  the  eigenvectors  (EOF  coefficients)  and 
compute  the  persistence  statistics,  would  take  at  most  a  few 
minutes,  and  could  actually  be  done  on  a  hand  held 
calculator  (Shaffer,  1982).  This  would  make  the  warning 
time  increments  equal  to  30,  '4  8  and  66  h  if  th^-  forecast 
technique  makes  use  of  the  analyzed  and  forecast 
geopotential  height  fields.  This  time  delay  problem  could 
be  partially  eliminated  by  using  the  previous  12-h  forecast 
fields  to  define  steering,  which  would  permit  the  technique 
to  be  available  at  the  time  JTWC  prepares  their  forecast  at 
00/12  GMT.  It  is  believed  that  little  effect  would  be 
apparent  for  the  shorter  forecast  periods,  although  there 
would  likely  be  a  larger  impact  on  the  extended  time  period 
forecasts. 
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II.        EOF    ANALYSIS     OF    SYNOPTIC    FORCING    APPROACH 
T0"'HECtIR71TURE 


This  research  is  a  sequel  to  a  synDptic  map  type  scheme 
of  Brown  (1981)  and  an  EOF  rap resentation  of  rhe  synoptic 
forcing  to  predict  tropical  cyclone  motion  in  the  western 
North  Pacific  Ocsan  by  Shaffar  (1932)  .  The  da-^a  fislds  and 
selection  of  storm  criteria  (see  descriptions  by  Brown, 
1981;  Shaffer,  1982;  and  Shaffer  and  Elsberry,  1932)  are 
exactly  the  same  in  all  three  studies.  A  short  summary  of 
the    fields   and   criteria    follow. 
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Figure   1.      Grid   Arrangement  [Jssi    in   Statistical 
-Climatological  Track    Prediction. 
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The  500  mb  D-values  are  sxtraczed  on  a  movable  120-point 
grid  (Fig.  1)  over  35°  in  latitjde  and  70°  in  longitude  with 
a  5°  increment.  The  tropical  ::yclone  is  always  centered  at 
grid  point  70  which  is  10**  from  the  sDQthernmost  boundary 
and    40°   from    the   western    boundary.  This    choice    is    made    -o 

include  both  a  portion  of  the  aid-latitude  flow 
characteristics  and  a  smaller  section  on  the  more  uniform 
tropical      height    fields      to  the      scuth   of      the    storm.  The 

movable  feature  is  incorporated  to  minimize  the  seasonal 
variation  that  would  occur  if  the  grid  was  stationary.  The 
earth's  curvature  (convergence  of  the  meridians)  is  not 
included  over  the  grid  but  oould  be  easily  added  for 
operational  use.  The  D-values  on  the  standard  63  x  63  grid 
from  uhe  FNOC  archives  are  deviations  from  the  st-andard 
atmosphere  geopotential  and  ara  fitted  to  the  120  point  grid 
by    a    Bessel    linear    interpolation    scheme. 

Reguireraents  for  •^rcpical  cyclones  to  be  included  in 
this    study    are: 

•  The  tropical  cyclone  must,  be  of  at  laasr  trooical  s-orm 
intensity  (>17  m/sec)  and  be  located  to  th.^  west  of 
180°; 

•  The  cyclone  must  exist  at  tropical  storm  intensity  for 
at    least    30    hours  as   analyzed   by   JTWC    (Guam) ; 

•  The  cyclone  must  be  located  between  10'^N  and  25°N  (This 
is  done  to  prevent  the  grid  from  antaring  the  southern 
hemisphere  and  to  avoid  possible  problems  arising  from 
the    convergence    of    the  meridians   at   higher    latitudes) ; 
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•  Ths  storms  must  exist  at  OD/12  GMT  for  the  scheme  to  be 
able   to   use   the    correspDnding   FNOC   data   analyses; 

•  2ach  storm  observation  included  mast  have  36  hours 
between  a  prior  cbservatioti  of  the  same  storm  to  insure 
pseudo-independence  within  the  sample.  This  very 
important  restriction  insures  that  the  statistical 
results   are   unbiased. 

One   approach  to   assess   spatial   variations    in   geophysical 

fields  is    through     the  use  of    Empirical      Orthogonal   Function 

(EOF)    analysis    (also     known   as    Principal   Component   analysis). 

EOF    analysis      may    be    applied      to    any      multi-dimensional    data 

set    to  separate   signal  from  noise,      and   identify   significant 

spatial      correlation    patterns       (Hotelling    (1933);         Kutzbach 

(1967);      Morrison     (1967);       Barnati:      (1978);       Richman     (1980); 

Shaffer    (1982)    and   others).        For   geophysical    use,      the    data 

are      placed      on      the      moveable   grid      with      each      grid      point 

representing   a   dimension    in  space.      The   iits.    are   represented 

as    an    .^xN    matrix,      with    M    the      number    cf    grid   points,       and    N 

the      number  of      data      cases.         Each      column      of    this      matrix 

represents  a    single    realization      of   the    observed    geophysical 

field.         An   MxM      matrix   R    is    taen   formed      by    multiplying    thB 

data    matrix   by      its   transpose.         The    form    of      R    depends    upon 

the      form    in      the    initial      data    matrix,         but    is      generally, 

either      a      covariance      or    correlation      matrix      of      the      grid 

points.         If   cne   wishes    to   give    equal    weighting   to   each    grid 
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point,  the  correlation  matrix  is  formed.  '  If,  on  the  other 
hand,  the  desire  is  to  give  greatest  weighting  to  those 
dimensions  (grid  points)  with  the  largest  univariate 
variation,    the  covariance    matrix    is   forms! . 

The  choice  of  correlation  versus  covariance  matrix  is 
critical  for  meteorological  fields  which  cover  a  large 
spatial  domain.  For  example,         in    analyzing      geopotential 

fields,  mid-latitude  grid  points  have  much  greater 
variability  than      tropical  grii    points.  If   the   aim      is   to 

focus  on  mid-latitude  features,  the  covariance  matrix  form 
of  R  is  preferable.  If,  on  tha  other  hani,  the  wish  is  -c 
de-cermine  variability  patterns  over  th=  entire  domain,  the 
correlation   matrix   form      of  R    is    favored.  The   cccuela-ion 

matrix  approach  mutes  the  EOF  response  in  higher  variability 
regions  (mid-latixudes) ,  whila  amplifying  the  response  in 
regions   of    lower   variability.  Both    approaches    were    tested 

in    this   study.  With  the   covariance    matrix      approach,      t-he 

total  explained  variation  in  the  D-valua  fields  for  the 
first  10  eigenvectors  was  more  than  94  76,  which  was  9% 
greater      than      the      correlation      eigenvsctors.  Regression 

analyses  from  the  covariance  method  were  shown  to  be 
inferior  to  those  from  the  correlation  msthod.  Thus,  tha 
correlation   matrix    method    is    usad   in    this    study. 
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After  forming  the  matrix  R,  minimization  of  the  total 
variability  leads  directly  to  the  eigenvalue  problem 
(Morrison,  1967).  The  eigenvector  associated  with  the 
largest  eigenvalue  is  the  single  vector,  out  of  all  possible 
vectors  in  space,  that  minimizes  the  summed  variation 
(squared)  of  all  of  the  observation  vectors  in  M  space. 
Furrhermore,    this  eigenmode  explains   exactly 


X, 


M 

E  X 

i  =  1 


of  the  total  variation  in  the  observations,  where  A  is  the 
eigenvalue  of  the  ith  mode.  Since  the  trace  of  a 
correlation  matrix  is  identically  equal  to  the  order  of  the 
matrix,  the  percent  of  total  variation  explained  by  any 
eigenvector  is  simply 

X, 


The  elements  of  each  eigenvector  may  be  loosely  interpreted 
as  the  multiple  correlation  of  the  corresponding  grid  point 
with  all  other  grid  points,  simultaneDusly.  The  exact 
meaning  of  the  eigenvector  elements  (loadings)  is  given  by 
Morrison     ( 1967)  . 
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The  aim  of  this  project  is  to  investigate  spatial 
variation  structures  of  analyzed  500  mb  D-value  fields  for 
the  western  North  Pacific  Ocean,  Therefore,  we  chDose  the 
correlation      matrix   approach.  The      initial      step   in      data 

processing   is   to   point-normalize    the    data.       Specifically,    we 
use 

N 

Dbar     (i)     =      J_      E      D(irj) 
N        j-1 

D*  (i,  j)      =      D  (i,j)    -   Dbar(i) 


[y^    (   D  (i>  J)    -    Dbar(i)     )**2      2 


J=1  N-1 

where  -j  is  the  map  (case)  index  and  i  is  the  grid  point 
index.  Thus,  D*  (i»j)  is  simply  a  QDrmally  distributed 
random  variable,  with  a  mean  at  each  grid  point  over  all 
maps      of       zero      and       a     variance      of      1.  This      allows      an 

equally- weighted,  grid        point      interpr = taticn        Df        the 

eigenvectors.  For      this      study,      the      lata      analysis      was 

completed    using   the   5  00    mb   fields   only. 

The      data    consist      of      D-values      of    ge opotential,  (the 

difference  between  the  observed  geopDtential  and  the 
standard        gecpotential,         so         that        a      positive        D-value 
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corresponds  to  higher  than  standard  gaopotentials)  ,  in 
meters.         EOF   analysis  was  performed    for    1 20-dimension    space 

(8x15      grid)       in      Fig.      1-  Examples      o£    rhe      first    9      EOF 

eigenvectors  (correlation  case)  are  shown.  in  Fig.  2.  It 
must  be  stressed  that  the  negative  of  the  pattern  shown  in 
these      figures   is     also   a      solution   to     tha    problem.  (For 

examples      of      700inb      and        850mb      eigenvectors      see      Shaffer 

(1982)  ,) 


Figure  2.   Eiaenvectors  1-9  Multiolied  by  130,   The  X^ 

Indicates  the  Lccarion'of  the  Tropical  Cycxone, 
Grid  Point  70. 

The   number  of   eigenmodes   to  retain   in   each  c3iS^      is 

determined   by  comparing   the   aigenvalues   of  the   analyzed 
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geopoteatial  fields  with  the  sigenvalues  of  a  Mon-e  Carlo 
simulation  of  random  fields  as  used  by  Preisendorf er  and 
Barnett  (1977).  The  eigenmodes  retained  are  those  that  are 
larger  than  the  random  eigenvalues  plus  3  standard 
deviations.  By  using  this  approach  we  are  99%  confident 
that  zhe  only  eignvectors  rstained  ars  those  containing 
actual  signal,  and  that  we  are  discarding  noise  only.  Using 
this  selection  procedure,  the  actual  nunber  of  significant 
eigenvectors  (EOF's)  retained  is  limited  to  10.  This  number 
of  EOFs  explains  approximatsl/  85%  of  ths  variation  in  the 
fields.  Table  I  outlines  the  increase  of  axplained  variance 
as  the  number  of  significant  eigenvectors  retained  is 
increased.  Fig.  3  is  the  graphical  reprBSsntation  of  Table 
I. 

The  total  sample  of  the  50U  st.orms  during  the  1967  to 
1976  time  period  is  divided  into  a  dependant  samplr  of  U5U 
cases  and  a  50-case  independent  sample  for  use  as  the  test 
sample.  This  is  the  identical  sampia  of  50  randomly 
selected  independent  storms  that  Shaffer  (1982)  tssusd  for 
displacement  forecasts.  Shaffer  (1982)  and  Shaffer  and 
Elsberry  (1982)  showed  that  for  rhis  independent  sample  of 
Tropical  cyclones,  the  forecast  error  froi  a  technique  based 
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TABLE    I 


Of 


Eigeavalues  and  Cumulativs   Pircent    Explained   Variance 
the  Normalized   D- value    Fields    at  500    mb.     (Shaf f sr , 1 932) 


Mode 

1 

2 

3 

H 

5 

6 

7 

8 

9 
10 
15 
20 
4  0 
60 
120 


Eigenvalue 
39.36 
21  .60 
9.29 
7.43 
6.08 
5.29 

u.oo 

3.13 
2,49 
2.11 
1  .03 

.58 

.11 

.04 

.00 


Explained  Variance 
33 
51 
59 
65 
70 
75 
78 
81 
83 
85 
91 
94 
98 
99 
100 


on  an  EOF  representation  of  -he  synoptic  forcing  is 
decreased  by  an  average  of  17%  <<  hen  compared  to  the  Dfficial 
JTWC  (Guam)  forecast  error. 
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Figure  3.       Eigenvalues  of    the    50  0    mb    D-value    Fields 
(circles)     and    the    ;ionte   ZarlD    j  =  r.era-ced 
Eigenvalues    (triaaglss)    Based    oa    100    Simulations 
of    Random    Data. 
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III.       DEFINITION    DF    RECURVArURE 

The  definition   of     a    recurvaturs    case   can      be   thought   of 

as    a   set   of   restrictions    based      on  the   obssrved   track   of   the 

tropical    cyclone.  The    definition      of   recurvature      that    is 

used   in   this      research  is    rhat   described      by   Leftwich    (1978, 

1982): 

"Recurvature  is  a  net  displacement  northward  of  315° 
during  the  forecast  period  or  the  attaining  of 
northeastward  morion  during  the  12  hours  prior  to  the  end 
of   the    forecast   period". 

Storms  are  net   considered    if   their   past    12    hour'motion    is    to 
the      east..  k      graphical   description      of       recurvature      and 

non-recurvaturs   cases    is    given   in   Fig.    1. 

Some  of  the  cases  present  rather  iifficult  forecast 
situations.  Figure  5  depicts  Typhoon  Trix  (initial  position 
at  OOGMT  21  August  1971),  which  for  tha  initial  36-h  period 
would  be  considered  to  have  recurved,  whereas  for  the  54- 
and  72-h  period  it  would  be  cDnsidered  as  a  non-recurvat ure 
case.  Figure  6  depicts  Typhoon  Polly  (initial  position  at 
12GMT  29  August,  197U),  which  for  the  36-h  period  is 
considered  to  have  recurved,  but  then  moves  westerly  to  a 
5U-h    position      which    is    considered      as   having      not   recurved. 
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Figure  ^.      Examples   of  Idealized    Eecurvature     (1,2),    Non- 
recurvature    (3)    and   a    Not    Considered   Case    (4) 

and      finally   moves      north-norrhw est   to      bs      considered    as      a 

recurvature  case  at   72  h.      Typhoon  Cora    (initial   position   of 

12GMT   02    October    1975)      shown    in    Fig..         7    is    an    example    of    a 

recurvature      case   a~      all   three      periods,       whereas        Typhoon 

Alice    (initial   position   at    12GMr    16    Septeubsr    1975)       in    Fig. 

8,    is  a   non-'recurvature   case    for    all   three    periods.      Typhoon 

Anita    (shown      in   Fig.      9)       would      net   bs    considered      for    the 

12GMT    22    July    to   OOGMl  23      July    initial   positions    due    to    the 

eastward      motion   of      the      tropical      cyclone,       but      would      be 

considered    for   the    12GMT    23  July    to    OOGIIT    24    July   period. 

A     summary   of      the   recurvature      statistics      for   the      SOU 

storm   sample      is   given      in  Tabla    II,  The    6      h    recurvature 
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Figure  5.      Typhoon   Trix   Case    Illustrating    an    Apparent 

Recurvatur€  Case  at    3  6   h,    and    Non-recurvature    at 
5U    and   72    h 
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Figure  5 . 


Tyohoon   Polly    Case   Illustrating    Pecurvarurs 
ana   72   h   and    Non-r sea rvature    at    54    h 
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Figure  7.   Typhoon  Cora  Case  Illustrating  a  Common 
Recurvature  Case 
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Figure  3.      Tvphoon    Alice    Case    II  lust  raT:ing    a    Non-racurvat  ur  3 
Case 

statistics      were  obtained      in      the      sams    nanner      as      defined 

earlier,    except   'hat    the    mcxion    during   the    past    6    h   »rfas    used 
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Figure  9,   Typhoon  Anita  Casa  Illustrating  the  Exclusion 
rhe  12GMT  22  July  zo    OOGMT  23  July  Periods  Du; 
the  EaSuward  Motion 
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for  The  second  consideration  fDr  recurvature  vice  the  past 
12  h  motion.  All  past  meridional  and  zonal  motions  were 
obtained  from  the  warning  positions  of  each  individual 
storm.  This      was      done      to      simulate      actual      operational 

conditions    (Shaffer,     1982). 

A  graph  of  the  recurvature  statistics  shown  in  Table  II 
is  presented  as  Fig.  10.  There  is  a  slow  increase  of 
probability  of  recurvature  with  time.  This  increase  may  be 
explained  by  the  following  reasoning.  A  tropical  cyclone  in 
its  developmental  stages  is  mora  often  in  the  easterly  trade 
wind   regime      and   moving    west      of    northwest       (non-recurver    as 
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TABLE    II 


Number  Df   Recurving    (rev).    Non-recurving    (nrcv)  ^    Nq   History 
(nohist)  ,    Not   Considered    (notcon)    and   tfie    Probaoility    of 
Recurvature    (prob)    for  the  504   Tropical   Cyclone   Sample 


tau 

rev 

nrcv 

nohist 

notcon 

prob 

06 

131 

234 

0 

139 

.359 

12 

125 

240 

0 

139 

.342 

18 

130 

235 

0 

139 

.356 

24 

134 

231 

0 

139 

.367 

30 

135 

230 

0 

139 

.370 

36 

158 

196 

22 

128 

.446 

42 

119 

169 

1  15 

101 

.413 

48 

125 

163 

1  15 

101 

.434 

54 

133 

154 

1  16 

101 

.463 

60 

133 

139 

1  42 

90 

.489 

66 

96 

1  11 

2  30 

67 

.464 

72 

103 

104 

230 

67 

.498 

76 

107 

100 

2  30 

67 

.517 

84 

99 

87 

2  54 

64 

.532 

defined)  .  As  the  tropical  cyclone  intensifies  and  moves  to 
higher  latitudes,  it  begins  to  interact  wirh  the 
mid-latirude  westerlies  and  therefore  is  more  likely  zo 
recurve.  Cheng-Lan  and  Saalar  (1982)  describe  the  racsi 
freguent  track  of  a  recurving  -ropical  cyclone  as  "=  regular 
clockwise  parabola".  That  is,  as  the  -^.ropical  cyclone  moves 
westward  it  also  moves  nor-hward ,  wirh  =  greater  likelihood 
of    being    north   of   a   3  15°    line    from   the    initial    position. 
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Figure   10.       Probability   of   Recur vatare    (504    Case    Sample 
Climatology) 

The   statistics      contained    la    Table      II    and   Fig,         10    are 

designated      as    -he      sample   clima tological      probabilities      of 

recurvature.         It      should    be    mentioned   thit      the    independent 

sample  statistics    for      recurvature   for   all    time      periods   are 

within  7%   of   the   sample   average    probability,    which   indicates 

consistency   between    the   two   samples.  These    values    will    be 

used   later   as   a   standard      for   evaluation    of    various   forecast 

schemes      in        predicting      recurvature        at      specified        time 

intervals.         That    is,    a   scheme   will   be   said    to   have   skill    if 

it        predicts        recurvature        more        accurately        than        the 

climatological   probability. 
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IV.       3EGRESSI3N    ANALYSIS 

Neumann      and        Lawrence      (19  75)         damonstrated        that      a 
St atistical- dynamical  modsl  using  prognostic  and 

persists nce-related  variables  in  a  aoQlinear  regression 
scheme  shewed  much  promise  in  forecasting  Atlantic  ocean 
tropical      cyclone      motion.  Shaffer    (1932)         also      used      a 

st atisuical-dynamical  scheme  that  used  10  20F  coefficients 
to  represent  the  geopctential  field  (analysis  time  only)  and 
a  linear  regression  scheme  to  predict  tropical  cyclone 
motion.  The        30F      analysis         not      only         simplifies      the 

representation  of  the  synoptic  forcing,  it  also  reduces 
significantly  the  cost  of  producing  an  accurate  and  viable 
forecast . 

The  importance  of  the  nonlinear  terms  in  regression 
schemes  has  been  demonstrated  by  Neumann  and  Lawrence  (1975) 
and  Paegle  and  Haslam  (1982).  although  noalinear  terms  were 
included  in  those  studies,  economic  considerations  and  the 
simplicity  of  application  of  a  linear  model  led  to  its 
adoption    in  this   study. 
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To  test  the  effect  of  iricrluding  EDF  coefficients  as 
predictors  of  tropical  cyclone  motion,  thres  separata  trials 
for  each  of  the  prescribed  times  were  completed.  The  first 
test  (refered  to  as  PERS)  was  accomplished  using  only  the 
persistence-related  parameters  and  the  initial  latitude  and 
longitude.  A  description  of  all  potential  predictors  is 
contained      in      Table      III.  (All      latitude      and      lor.gi-uie 

displacements  (displ)  are  in  nautical  miles.)  The  second 
experiment         (referred     to        as      SOFS)  utilized      10        EOF 

coefficients  and  the  initial  latitude  and  longitude.  The 
third  experiment  used  a  combination  of  the  persistence  and 
EOF  predictors  ar.d  is  noted  as  P&E.  It  should  be  mentioned 
that  the  cyclone  maximum  sustained  wind  speeds  at  12,  24  and 
36  h  prior  tc  the  initial  observation  -4=rz  also  incorporatad 
as    possible   predictors.  This    was   done    in      accordance    with 

Hiehl's  hypothesis  (197  1)  of  maximum  intensity  being 
associated  with  recurvature  for  Atlantic  tropical  cyclones. 
However,  none  of  these  additional  predictors  were  chosen 
during  the   stepwise   regression. 

The  504  case  tropical  cyclone  sampla  was  analyzed  for 
recurvature  as  described  in  Chapter  III.  A  value  of  1.0  was 
assigned    to   the    storms  that  reourvei      and    a    value   of    0.0    was 
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TABLE   III 

Potantial   Predictors    Used   to    Develop   the    Regression 
Equations   for    Each   of  the    Three    rime    Intervals 

Potential    Predictor        Name  Desgr j.pt ion 

Variable  Number  (Coef icient-Coef ) 

1  C0F1  Eigenvector     1 

Coefficient 

2  C0F2  Eigenvector     2 

Coefficient 

3  C0F3  Eigenvector     3 


Eigenvector 
Coefficient 


a  COFU  Eigenvector    U 

Coefficient 

5  C0F5  Eigenvector    5 

Coefficient 

6  C0F6  Eigenvector    6 

Coefficient 

7  C0F7  Eigenvector    7 

Coefficient 

8  C0F8  Sigenvecror     8 

Coefficient 

9  C0F9  Eigenvector    9 

Coefficient 
10  COF10         Eigenvector    10 

Coefficient 
11'  PLATO         Initial   Cyclone 

Latitude 

12  PLAT1  Past  12  hour 

meridional  displ 

13  PLAT2  Past  2U  hour 

meridional  displ 

14  PLAT3  Past  36  hour 

meridional  dis^l 

15  PLATU         Meridional  displ 

for  prior  12-2U  hrs 

16  PLATS         Meridional  displ 

for    prior    12-36    hrs 

17  PLAT6         Meridional  displ 

for  prior  24-36  hrs 

18  PLONO  Initial  cyclone 

longitude 

19  PL0N1  Past  12  hour 

zonal  displ 

20  PL0N2  Past  24  hour 

zonal    displ 

21  PL0N3  Past    36    hour 

zonal    displ 

22  PLONU  Zonal    displ 

for    prior    12-24   hirs 

23  PL0N5  Zonal  displ 

for  orior  12-36  hrs 

24  PL0N6  '    Zonal  displ 

for    prior    24-36   hrs 

assigned  to  the  cyclones  -hat  had  not  recurved  at  the  time 
of  interest.  A  value  of  999.  was  assigned  to  cyclones  that 
had    no  history   at    the   time   period   of    interest. 
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The  BMDP  stepwise  linear  rsgression  program  screened  the 
dependent  data  set  for  storms  with  past  12  h  longitude 
displacement  (PL0N1)  that  was  to  the  east  and  also  tropical 
storms  that  had  no  history.  This  screening  accounts  for  the 
different  number  of  cases  considered  for  the  regression 
routine  and  the  cliiratology  cases  at  the  specific  time 
periods  of  interest  (see  Table  II).  The  equations  at  36,  54 
and  72  h  (see  Tables  IV,  V,  7  1  below)  are  d-rrived  from  a 
stepwise  linear  regression  using  an  F  -  ratio  of  U.O.  The 
F-ratio  is  a  test  of  the  significance  of  the  independent 
variable  coefficients  in  the  regression  equations.  The 
number  of  cases  of  the  dependant  sample  that  are  selected 
for  each  regression  period  is  noted  in  the  corresponding 
table.  The  number  of  cases  dei^reases  with  time  due  to  the 
restrictions  for  inclusion  discussed  in  Chapter  II.  Of  note 
is  that  most  of  the  cases  at  72  h  are  a  subset  of  the  5U  h 
data  set,  which  is  in  turn  a  subset  of  the  36  h  data  set. 
The  R  Squared  Statistic  ,RSQ,  (percent  Df  variance  in  the 
predictar.d  explained  by  the  ragression  equation)  will  be 
used  for  discussion  of  the  results. 
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The  linear  regression  equations  produced  from  the 
potential  predictors  of  recurvatur=  or  non-recurvature  are 
of    the  form: 

?(0:1)    =   aO    +    alXI    +   a2X2    +   a3X3    + +    a7X7. 

Tables   IV,    V,   and   VI    list  the  coefficients    and   variables 

for    each      of   the      three   tests  generated    froi      the   dependent 

data     sets.        The     number      of  cases      used      to   generate        the 

coefficients   are   259,    206    and  13  8   respectively. 

TABLE    IV 

Regression   Equations   for    36    h   Recurvatur=    for    Each   of    the 
Three  Methods  Tested    (Dependent    Sample) ,    and    Mean  Values 
and   the    Standard    Deviation  from    the    Mean     (STD) 


PEES 

SOFS 

PSE 

MEAN 

STD 

intercept 

.  1852 

^.2871 

.5502 

PLATO 

.0167 

19.4 

5.2(0) 

PLONO 

.0054 

133.2 

12.9  (O) 

PLAT2 

.0017 

.0014 

96.4 

93. 2  (NM) 

PL0N1 

.0031 

.0031 

-92.6 

55.3  (NM) 

COF     1 

-.0326 

-.0215 

.81 

5.  12 

COF    3 

.0367 

-.24 

2.89 

COF    5 

-.0302 

.20 

2.34 

COF    7 

.0341 

,02  99 

-.07 

1.91 

COF    8 

.0493 

-.20 

1.77 

COF     10 

-.0545 

-.03  97 

-.  10 

1.48 

Each    regression   coefficient      is    unique   excap*    for      the    value 
of    PL0N1    for   the   PERS    and    P8S      test   at   36    hours.         It    should 
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TABLE 

V 

Regression 

Equations 

for  54  h  1 
Table 

R ecurvaturs 
IV 

as  Descr: 

Lbed  in 

PERS 

EOFS 

PSE 

MEAN 

STD 

INTERCEPT   . 

.uuna 

.4957 

.4873 

PLAT2 

.0019 

.0014 

95.6 

86.6 

.0022 

.0016 

-9  3.0 

56.6 

COF  1 

-.0312 

-.0252 

1.23 

4.8 

COF  3 

.0371 

.0237 

-.53 

2.74 

COF  8 

.0618 

.0465 

-.21 

1.82 

COF  10 

-.0514 

.05 

1.50 

TABLE 

VI 

Regression 

Equat 

ions 

for  72  h  1 
Table 

R ecurvaturs  as  Descr 
IV 

iDed  in 

PZRS 

EOFS 

P&E 

MEAN 

SID 

intercept 

.5487 

.5311 

.5821 

PLAT4 

.0028 

.0026 

44.  1 

48.  1 

PL0N1 

.0027 

.0024 

-89.  1 

52.2 

COF  1 

-.0  3  24 

-.0230 

2.02 

4.55 

COF  3 

.0410 

-.75 

2.77 

be  noticed  that  the  value  of  PL  ONI  must  necessarily  be  less 
than  zero,  because  a  PL0N1  greater  than  zero  represents 
easterly  movement,  which  is  considered  as  having  already 
recurved . 

The  total  number  of  combinations  of  predictors  and 
regression  coefficients  will  not  be  discussed,  but  several 
combinations    will   be    mentioned.         PLAT2      (past    24    h   latitude 
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displacement)  is  almost  always  greatar  than  zero  (the 
exception  being  trochoidal  oscillating  stDrras) ,  and  makes  a 
considerable  contribution  towards  recurvature  (pradictand 
close  to  1.0)  especially  for  large  displacements-  PLATO  and 
PLONO  (initial  latitude  and  longitude)  are  always  greater 
than      zero,  with      larger       values      being        favorable      for 

recurvature      prediction.  The      coefficient      of      the      first 

eigenvector  (COF  1)  has  the  opposit.e  effect.  Physically, 
the  stronger  the  zonal  pressure  gradient  (large  COF  1),  the 
lower  the  probability  of  recurvature.  The  opposite  is  also 
true;  the      probability      of        recurvature      is      higher      for 

extremely  weak  zonal  pressure  gradients  as  illustrated  in 
Fig.    11. 

Analysis  of  the  sample  mean  values  of  the  coefficients 
of  EOF  1  (assuming  a  normal  i istr ibution  about  the  mean) 
indicates  that  the  number  of  recurving  tropical  cyclones 
would  be  less  than  the  number  of  non-recurving  storms.  For 
example,  the  number  of  recurving  storms  at  36  h  (Table  II) 
is  158  and  the  number  of  non-r acurv ing  storms  is  196.  The 
coefficient  of  EOF  5  is  also  quite  interesting.  The  more 
the  EOF  5  pattern  describes  the  synoptic  forcing  (the  larger 
the    value      of    COF    5)  ,         the   mora      likely    the    storm      will    not 
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Figure   11. 


750 


X 

fe    500 
o 


250 


X        0 

I- 

8  -250 


:COF  1  BETWEEN  -9  AND  -30 
I  I  I  I  I  I  I  I  I  I'  I  I  I  I  I  I  I  I  I  I  I  I 


-250     0        250     500     750 
WEST  EAST 


q: 
o 


500  - 


250 


-    b 


8-250 


-COF  1  BETWEEN +30  AND +9 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  t  t  I  I  I  I  I  I 


-750  -500  -250     0        250 
WEST  EAST 

storm   Tracks    Relat,ive   to    the    I'nitial   Position 
for    All   Storms   with    an  Orthogonal    Coefficient    1 
that    is:       (a)     Bet  wean   -9    and   -33;    or    (b)     3e-ween 
+  9   and   •••3  0.      Each   Cross  Indicates   a    l2-ii 
Displacement     (Sha  ff  er,  1982) 


recurve.  The  converse  is  also  true  in  that  the  stronger  the 
inverse  or  the  negative  of  the  pattern  (a  strong  shortwave 
between  two  ridges) ,  the  more  likely  the  cyclone  is  -o 
recurve. 

This  is  by  no  means  a  complete  discussion  of  all  the 
possible  54  linear  combinations  of  predictors,  coefficients 
and    their    inverses.        It    does      show  how   a    simple   combina-ion 


37 


of  EOF  patterns  can  be  applied  with  straight-forward 
meteorological   reasoning. 

A.       36    H    REGRESSION    CASE    ANALYSIS 

As  expected,  the  predictors  selected  for  the  PE3S 
regression  analysis  included  the  past  12-h  longitude 
displacement,  the  past  24-h  latitude  displacement  and  the 
initial  latitude.  Simple  extrapolation  techniques  do  well 
for   these      relatively   short-tiaie    forecasts.  The  longitude 

displacement  over  the  past  12  hours  reduces  the  variance  in 
the  predictand  (increases  RSQ)  slightly  more  than  the 
latitude  displacement.  The  selection  of  initial  latitude  as 
a  predictor  is  consistent  with  the  variation  of  the 
probability  of  recurvature  with  time  (and  therefore  latitude 
for  an  idealized  storm  track)  as  discussed  in  Chapter  II. 
The  total  explained  variance  of  the  predictand  for  the 
persistence-related   predictors   <*  as  29.2?5. 

The  individual  EOF  predictors  reduce  less  of  the 
variancs  than  the  persistence-type  predictors  (Table  VII). 
Notice  that  eigenvector  1  (which  is  an  indicator  of  the 
strength  of  the  zonal  pressure  gradient  and  large  global 
wave  number)  causes  the  largest  reduction  in  the  variance  of 
the    predictand.         The   remaining      five    eigenvectors   that    were 
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selected  are  indicative  of  smaller  scale  forcing  features, 
and  reduce  the  variance  of  the  predictand  much  less  that  the 
first  eigenvector.  Shaffer  (198  2)  indicated  the  possibility 
of  weak  northwest  motion  (or  southeast  motion  for  the 
inverse  mode)  would  fce  associated  with  eigenvector  7  and  a 
strong  meriodional  forcing  being  associated  with  eigenvector 
3.         Of      special   interest   is      eigenvector    5.  Al-hough   the 

reduction  in  variance  associated  with  this  eigenvector  is 
small,  its  inverse  pattern  appears  similar  to  a  strong  short 
wave  in  the  westerlies  which  could  advect  the  srcrm 
northward.  The     last      predictor     selected,        the      initial 

longitude,  again  is  consistent  with  the  variation  of  the 
recurvature      probability      with      time.  The      total      of      the 

explained  variance  for  the  S0?5  test,  was  26.5?b  ,  which  is 
2.7%   less   than   using    persistence    only    predictors. 

The  selection  of  predictors  for  the  final  test  shows  a 
combination  of  both  P ERS  and  ED FS  predictors.  The  largest 
reduction  of  the  variance  occurs  with  past  motion 
(persistence  related)  parameters.  The  strength  of  the  zonal 
pressure  gradient  appears  as  the  next  most  important 
predictor  with  smaller  seals  features  (EOFS  10  and  7) 
contributing   small    increases   to    the   explained   variance.      The 


39 


total  of  the  explained  variance  of  the  prsdictand  is  3a. 0% 
which  is  greater  than  either  the  PERS  or  EOFS  methods,  but 
still  only   a   small   fraction  of  -he  total   variation. 

Part  of  the  reason  thaz  more  of  the  predidictand 
variation  could  not  be  explained  may  be  related  to  the 
arbitrary  definition  of  recurvature,  e.g.  a  movement  10 
320^    is    recurvature    (value    1.0)  whereas    310°    is   not     (value 

of  0.0).  When  a  yes/no  decision  is  applied  across  a  line  as 
in  this  situation,  one  cannot  expect  an  explained  variance 
near  100%.  The  methcd  could  still  provide  useful  guidance 
for  the  forecaster  if  the  percentage  of  correct  forecasts  is 
higher  than  other   methods   available. 

The  entry  of  PLATO,  and  then  the  subsequent  deletion  as 
the  EOF  information  is  added,  is  in  indication  cf 
overlapping  information.  Notice  (Table  VII)  that  PLATO  dees 
not    enter    in   scheme    (c)    at  all. 


40 


TABLE    (^11 


Summary   of    Stepwise   Linear   Regras^ion    Results    for    36    h    (259 
cases)    R    Squared    (RSQ),    Increase    m   R    Squared    (RSQ   Incr) 
and  ?   to   Enter  or  R  emove 


a) 

PEESistence 

only 

m  ^  ^m   ^  ^B    ^m^m    ^   ^  ^b 

Step 

)          Parameter 

ESQ 

RSQ    Incr 

F    enter 

F   remove 

1 

2 

3 

PL0N1 

PLAT2 
PLATO 

.  151 
.263 

.  292 

.151 
.1  12 
.0  29 

45.65 
39.06 
10.42 

b) 

FOFS    only 

1 
2 
3 

4 
5 
6 

7 
8 
9 

C0F1 

PLATO 

C0F8 

C0F7 

C0F3 

COF10 

PLATO 

C0F5 

PLONO 

.  121 
.  161 
.  187 
.  206 
.  227 
.  24  3 
.  234 
.  248 
.  265 

.1  21 

.0  40 
.0  26 
.0  19 
.021 
.0  16 
-.009 
.0  14 
.017 

35.40 
12.04 
8.20 
6.04 
6.85 
5.47 

4.54 
5.75 

2.94 

c) 

Persistence 

& 

SOFS 

1 

2 
3 
4 
5 

PL0N1 

PLAT2 

C0F1 

CCF10 

C0F7 

.  151 
.  263 
.  313 

.  327 

.  340 

.151 
.1  12 
.050 
.0  14 
.0  13 

45.65 

39.06 

18.54 

5.  10 

5.  18 

B.        54    H    REGRESSION    CASE    ANALYSIS 

The  selection  of  PERS  predictors  at  54  h  (Table  VII)  are 
the  same  as  at  36  a  (Table  71)  except  for  the  dele-ion  of 
the  initial  latitude.  Had  th9  F-rest  been  set  for  3.9  vice 
4.0,  the  initial  latitude  would  have  been  selected 
subsequent  to  PL0N1,  as  was  th=  case  for  36  h.  At  54  h  the 
latitudinal  displacement  for  the  past  2^  h  decreassd  the 
variance  of  the  predictand  niore  than  iid  the  past  12  h 
longitude    predictor.       These  two    predictors    comprise    16.4%   cf 


41 


the  variance  of  the  predictand,  which  is  almost  half  the 
value  obtained  at  36  h.  This  decrease  is  not  surprising, 
because  simple  extrapolation  techniques  have  less  validity 
at    5a   h   than   at    36    h. 

The  number  of  predictors  selected  in  the  EOFS  test 
decreased  to  four,  compared  to  seven  at.  36  h.  These  four 
predictors  at  54  h  are  all  included  within  the  first  five 
predictors  a-  36  h.  The  north/souuh  pressure  gradient  (EOF 
1)  '  again  contributed  to  the  largest  reduction  in  the 
variance  of  the  predictand  (see  Figure  2).  EOF  8  appears  to 
be  similar  to  EOF  5  in  that  it  is  associated  with  a 
mid-laxi tudfe  short  wave  combined  with  an  area  of  high 
pressure  directly  south  of  rhe  tropical  storm.  EOF  8  is 
correlated  directly  with  recurvature  as  are  EOFS  1  and  3. 
EOF  10  is  indirectly  correlated  with  recarvature  due  to  the 
high  pressure  area  to  the  northwest  of  the  storm  center. 
The  stronger  this  pattern,  the  less  likely  it  is  that  the 
storm  will  recurve.  The  total  of  the  explained  variance  of 
the  predictand  is  IS.Ufi  which  is  smaller  than  the  value  for 
36  h  using  EOFS  only,  but  it  is  larger  than  54  h  regression 
using   persistence      parameters   only.  This   result      suggests 

that   synoptic    parameters    begin   to   explain    more   of   the    motion 
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at  later  times  than  does  the  extrapolation  associated  with 
persis-ence   variables. 

The   combination   of     the  two    types   of      predictors   explain 
24.6??     of    the      variance    of     the    predictand.  The  only      EOF 

predic-or  that  does  not  appear  in  P&E  is  eigenvector  10. 
The  inclusion  of  the  eigenvectors  as  predictors  increased 
the    RSQ   amount  by   8.2^.  The   decrease    from    34.0%   explained 

variance  at  36  h  to  24.6%  at  54  h  was  expected  due  -o  the 
likelihood  of  more  tropical  cyclone  interactions  with  the 
synoptic    environment. 

TABLE   VIII 

Summary    of    Stepwise   Linear    Ragression    Results    for   54    h 
(206    cases)    as    Descrioed  in   Table    VII. 

a)  PEHSistence   only 

Step  Parameter  ESQ  RSQ    Incr  F    enter        F    remove 

1  PLAT2  .102  -1  02  23.  26 

2  PL0N1  .  164  .062  15.06 

b)  EO?S    only 

1  C0F1  .  C74  .074  16.20 

2  C0F8  .  127  .053  12.42 

3  C0F3  . 170  .043  10-41 

4  COF10  .194  .024  6.03 

c)  Persistence   &   SOFS 

1  PLAT2  ,102  .1  02  23.26 

2  PL0N1  .  164  .062  15.06 

3  C0F1  .203  .039  9.93 

4  COFB  .230  .  0  27  6  .  97 

5  C0F3  .246  .0  16  4.15 
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C.        72    H    REGRESSION    CASS    ANALYSIS 

The  selection  of  predictors  for  ths  72  h  prediction 
includes  a  new  predictor,  PLAT  4  (latitude  displacement  for 
the  period  12  to  24  h  prior  to  the  initial  time).  This 
predic-or  is  inherently  positive  and  therefore  adds  to  the 
probability     of      recurvature.  The      past       12      h      longitude 

displacement  is  the  remaining  predictor  chosen.  The  total 
of  the  explained  variance  using  these  two  predictors  is 
15.05?.  Although  this  is  less  than  half  of  the  value  of  the 
36  h  regression,  it  is  only  slightly  lass  than  the  16.5% 
obtained    for   the   54   h   regression  . 


TABLE    IX 

Summary    of    St 
(138 

epwise   Linear 
cases)    as    Desc 

Regression    Results 
:ri5ed   in   Table    VII. 

for    72    h 

a)       PERSistence 

on 

ly 

Step          Paramet 

er 

BSQ 

RSQ    Ir 

icr 

F    enter 

F    removs 

1  PL0N1 

2  PLAT4 

.  074 
.  150 

.0  74 
.0  76 

10.84 
12.01 

b)       ZOFS    only 

1  C0F1 

2  C0F3 

.  076 

.  123 

.0  76 

.0  52 

11.27 

8.03 

c)       Persistence 

& 

SOFS 

1  C0F1 

2  FL0N1 

3  PLAT4 

.  077 
.  132 
.  192 

.0  77 
.055 
.  060 

11  .27 

8.63 
9.94 

The    EOFS-only      regression    produced   resul-s      that    explain 
less      of      the      variance      of      the      prsdictind      than      iid      the 
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persists nce-type  predictors.  Coefficients  1  and  3  explain 
12.8^  of  the  total  variance  of  the  predictand,  which  is 
almost  identical  to  the  value  for  the  first  steps  of  the  54 
h  regression.  However,  two  diffsreni.  coefficients  entered 
the    54  h    regression   eguation. 

The  combination  cf  two  groups  of  predictors  is  again 
superior      to      either        individaal     group.  At      72        h      the 

north/south  pressure  gradient  is  more  important  than  the 
longitudinal  displacements  84  h  earlier  (PL0N1)  or  the 
latitude  displacement  34  to  9  6  h  prior  t;o  the  initial 
position  (PLAT4) .  Use  of  extrapolation-ty ps  methods  can  not 
be  expected  to  contritute  as  much  as  the  synoptic  parameters 
at  such  long  time  intervals.  As  the  tropical  cyclone  moves 
into  the  mid-latitudes,  persistence  is  comparatively  less 
useful  as    a   predictor    than   in    the   tropics. 

The  amount  of  explained  variance  cf  the  recurvature 
event  decreased  with  time  for  each  of  the  three  test  groups. 
The  superior  method  (as  far  as  explaining  the  variability  of 
the  predictand)  is  the  combination  of  persistence-related 
variables  and  the  20F  coefficients  which  describe  the 
synoptic    patterns.  Although    the   PS S    technique   consistently 

explained    more   of      the   variance    than    the      other    two   methods. 
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further  nesting  is  necessary  to  determine  if  the  EOF 
predictions  are  superior  to  other  techniques  such  as 
climatology   or    mere    chance. 
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V.       TESTIS    THE   3S3RESSI0N    ANALYSIS 

Choosing  a  forecast  verification  scheme  is  very 
important  for  testing  the  ability  of  a  regression  model  or  a 
prognostic      method      tc      forecast      a      specific      event.  The 

geophysics  community  is  just  one  of  the  many  scientific 
communities  xhat  desire  an  accurate  and  acceptable  method 
for  describing  the  skill  of  a  given  model.  Muller  (19U4) 
provided  a  literature  survey  of  55  short-range  weather 
verification  methods  that  were  developed  during  1884  -o  1943 
in  six  different  countries.  Muller  divided  rhe  schemes  into 
two  areas:  comparisons  of  the  forecasts  and  observations  as 
a  ratio  (percentage  correct)  and  verification  methods  which 
are  compared  to  chance  or  climatology,  the  so-called  'blind' 
forecasts. 

In  the  early  verification  studies,  the  schemes  for 
comparing  forecasts  with  observations  had  to  be  simple  for 
the  ease  of  application.  Verification  :nettiods  today  compare 
the  prediction  skill  against  some  climatological  value  or  to 
pure    chance   as   a    point   of    reference- 
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Thrse  skill  scores  will  be  iiscussei:  the  Heidke  Skill 
Score  (HSS),  the  Brier  P-Scora  (BPS)  ,  and  the  percent;  of 
correct  forecasts  (PCF)  .  The  Heidke  skill  score  can  be 
written   as: 

HSS   =    7P    -C 


r   -  c 


where  7F  is  the  number  of  correct  or  verifying  forecas-s  of 
an    event    (ei-her   correctly     fcrscast    recurvers    (C?R)  ,      i.e., 

predictand  is  greater  than  or  equal  to  an  assigned 
climatology  probability  (Table  1 1)  ;  or  a  correctly  forecast, 
non-recurver  (CFNR)  ,  i.e.,  ths  predictand  is  less  than  the 
climatology  probability) ,  T  is  the  total  number  of  forecasrs 
made  and  C  is  aqual  tc  the  number  of  correct  forecasts  ihat 
would  be  expecxed  from  chance.  In  rhis  tsst  C  is  always 
one-half    of   the      total  number    oi    forecasts.  The   HSS    would 

have  a  maximum  value  of  1.0  for  perfect  forecasts  and  a 
value  of  zero  when  the  number  of  correct  forecasts  is  equal 
-o  the  number  of  forecasts  expscted  to  be  correci  by  using 
chance.  Negative  values  occur  if  the  method  is  inferior  to 
the    comparative    methods. 
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The  Brier  or  P-Score  has  been  widely  accepted  as  a 
method  to  test  probability  forecasts.  The  binary  form  of 
the    EPS   can    be    written  as: 

N 


S    =  J_   *    V^(FCST(i)     -    ACr(i))**2 


BP 

N  M 

where  FCST  is  the  regression  equation  pradictand  and  ACT  is 
either  1.0  if  the  storm  acta  ally  recarved  or  zero  if 
recurvature  did  not  occur.  A  perfect  3?S  is  a  value  of  zero 
and  the  worst  possible  forecast  weald  be  indicated  as  2.0. 
It  is  readily  apparent  that  the  BPS  is  actually  the  mean 
square  error   of   the  s  €t   of  probabilistic   forecasts. 

The   percent   of      correctly   forecast   avants    (PCF)         can    be 
written   as: 


PCF    =   CFH    +    CFNR 


where   CFR,      CFNR      and    T   are  the    same      variables    as    described 
earlier. 

As  discussed  in  Chapter  IV,  the  predicrand  values  should 
lie  en  -he  continuous  interval  between  1.C  and  0.0,  The 
unfiltered  predictand  statistics  for  both  the  dependent  and 
independent      samples    for      all    periods     appear      in   Table      XI. 
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Values  larger  than  1-0  (indicative  of  a  strong  recurvaturs 
signal)  and  smaller  than  0.0  (indicating  strong 
non-'recurvature  signal)  exist  as  maximum  and  minimum  values. 
These  values  outside  of  the  defined  range  are  due  to  the 
assumed  linear  relationship  between  the  predictors.  Table 
XII  provides  similar  statistics  for  the  predictands  after 
values  that  were  mor=  than  1.0  and  less  than  2.0  were  set 
equal  to  1-0  and  negative  values  were  set  equal  to  0.0, 
Comparison  of  the  dependent  and  independent  samples  of  bc-h 
tables  indicate  that  the  two  sets  are  quite  similar.  The 
different  number  of  cases  included  in  several  tests  is  due 
to  the  predictors  PLAT2  and  PL0N1  being  excessively  large, 
(purposely  assigned  flagged  values)  which  causes  several 
predictand  values  for  PERS  and  P&E  to  be  in  excess  of  2.0, 
the  upper  filter  limit- 
Tables  XI  and  X n  indicate  rhat  the  variance  of  the 
dependent  predictand  is  greater  than  or  equal  -o  the 
variance  of  the  independent  predictands  in  15  of  18 
comparisons.  This  is  evidently  due  t.o  tha  small  size  of  zhe 
independent  sample. 

It   also  should   te  noriced   that   the  dependent   sample 
variances  are  largest  at  36  h  and  decreasi  in  time.    In  all 
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dependent  data  set  comparisons,  the  variaaca  of  the  P&E  test 
is  larger  than  the  other  two  tests.  This  is  most  likely  due 
to  the  'signal*  of  either  recurvature  or  non-recurvature 
being  strongest  for  this  set  which  includes  both  persistence 
and  synoptic  features.  The  smaller  values  of  the  variance 
in  Table  XII^  as  compared  to  Table  XI,  is  due  to  the 
filtering   of   the  predictands,    which   is  described   above. 

TABLE    XI 

Predictand   Statistics    for   the   36,    54    and    72    h   Before 
Filtering.      The   Number    of  Cases,    Mean,    Maximum    and   Minimum 
Values    {max.min)    and    the    Variance    (var)    of    -he   Predicted 
Values   tor   "the   Dependent    and   Independent   Data    Sets   are    Shown 


Dependent    Set 
PESS         ECFS         PSE 


Independent    Set 
PEP.S         SOFS  PSE 


cases 

(3^9) 

^4^ 

QW 

36    h 

mean 

.49  8 

.432 

.512 

36    h 

max 

1.790 

1.  550 

1.810 

36    h 

min 

-.256 

-.  243 

-.2  59 

36    h 

var 

.120 

.  083 

.151 

cases 

(271) 

(345) 

(271) 

54    h 

mean 

.4  86 

.434 

.4  85 

54    h 

ma  X 

1.950 

1.  390 

1.780 

54    h 

min 

-.021 

-.151 

-.112 

54    h 

var 

.065 

,  G49 

.0  90 

cases 

<PJ) 

(235) 

^?J^^ 

72    h 

mean 

.485 

.  427 

.482 

72    h 

max 

1.36  0 

1.  06  0 

1.340 

72    h 

min 

-.037 

.  099 

-.038 

72    h 

var 

.053 

.  027 

.065 

(48) 

.45U 

1.  020 

-.259 

.096 

(49) 

.436 

1.  290 

.025 
.063 

{U8) 

.473 

1.050 

-.191 

.084 

.43  8 

.  867 

-.110 

.047 

.983 

-.082 

.041 

(^3) 
.45^ 

.889 

-.392 

.073 

(38) 

.459 

.  978 

-.  048 

.053 

(38) 
.4  34 
,821 
.121 
.028 

(33). 
.512 

1.070 

-.132 

.062 

The      summaries      of      the   HSS      which      compares      the 


-•n^Ao, 


methods  of  testing  using  the  climatology  probabilities  of 
recurvature  for  both  the  dependent  and  iniapendent  data  se-s 
are    presented    in   Table  XIII.         I-   should    be    no-iced    uhan    zhe 
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TABLE  XII 

Predictand  Statistics  fop  36,  54  and  72  h  After  Filtering  as 

Described  in  Table  XI. 


Dependent  Set 


Independent  Set 


cases 

PERS 
(349) 

ECFS 
(433) 

PSE 
(350) 

36    h 
36    h 
36    h 
36    h 

mean 
ma  X 
min 
var 

.486 

1.000 

.000 

.088 

.427 

1.000 

.  000 

.  C6  9 

.488 

1.0  00 

.000 

.099 

cases 
54    h 
54    h 
54    h 
54    h 

mean 
max 
min 
var 

(271) 

.47  8 
1.000 

,00  0 
.051 

1.  coo 

.  000 
.  046 

(271) 

.4  72 
1.0  00 

.000 
.068 

cases 
72    h 
72    h 
72    h 
72    h 

mean 
max 

min 
var 

(176) 

.483 

1.000 

.000 

.04  9 

(235) 
.427' 
1.  000 

.  C99 
.  027 

(176) 

.4  77 

1.000 

.0  00 

.057 

PERS 
(48) 


EOFS 
(49) 


P&E 
(48) 


.466  .430  .479 

1.000  1.000  1.000 

.000  .025  .000 

.082  .054  .075 

(431  ^^3)  (43L 

.443  .452  .467 

.867  .983  .889 

.000  .000  .000 

.042  .039  .058 


(38), 

.460 

(38) 

.4  34 

mi 

.978 

.821 

1.000 

.000 

.121 

.000 

.052 

.028 

.057 

number   of  cases  in   each  test   does   not   change  from   the 
earlier  values. 

TABLE  XIII 

Heidke    Skill  Score  Osina    -h=    Sample   Clima-::oloay   as    a 
Standard    cf  Comparison    for    Defining   aecurvature   Versus 
Non-recurvatura    for   'he   Fiiterad    Prsdicrands    (as   in    Table    II) 
A    HSS    Of   1.0   is    Perfect. 


Dependent   Set 

PERS         EOFS  P&E 

36    h  .553         .517         .526 

54    h  .505         .422         .565 

72    h  .431  .254         .454 


Independent    Set 
PERS         EOFS  P&E 

.464         .296  .469 

.442         .536  .440 

.431         .363         .510 


The  information  pertaining  to  the  skill  of  the  HSS  in 
Table  XIII  must  be  carefully  interpreted.  When  using  the 
HSS,      any      predictand   that      is   greater      t^an    or      equal    to      a 
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specific  value  (in  this  case  a  climatological  percantage  of 
recurvature  for  that  forecast  iaterval)  is  given  a  forecast 
recurvature  value  of  1.0.  Any  predictani  that  is  less  than 
the  climatological  value  is  considered  to  be  a 
non-recurvature  case,  and  is  assigned  a  value  of  0.0.  These 
•newly  assigned*  predictands  are  then  compared  to  the  act.ual 
storm   recurvature  values.  If    the   forecast    and     the    actual 

value  match,  a  verifying  forecast  is  tallied  (C?R  or  CFNR  is 
incremented  by  one). 

A  graphical  coiparison  of  HSS  versus  a  varying 
climatological  value  cf  probability  of  recurvature  (used  for 
comparison)  is  shown  in  Fig.  12.  The  maximum  values  in  the 
HSS  are  close  to,  although  nor  necessarily  associa-ed  with, 
the  climatological  probability  values  indicated  for  the 
three   different    time    periods. 

Another  point  needs  to  be  Mentioned  concerning  the  HSS. 
The  climatological  prcbabili-y  values  for  recurvature  (Table 
II)  at  the  specified  times  should  be  compared  to  -he  means 
and  variances  cf  the  predictands  shown  in  Table  XII.  All 
three  climarclogical  values  of  recurvatura  fall  within  the 
first  standard  deviation  around  the  mean  in  ail  xhree  tests. 
With   this    in   mind,      ore    would      predict,   thai    ae-hods   with   the 
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LEGEND 
0  -  35  H  PIE 
o  -  54  H  PiE 
A  -  72  H  PiE 


O.0SO.12S  0.200  0.275  0.350  0.12S  O.SOO  0.575  0.650  0.725  0.800  0.875  0.950 

PROBABILITY  OF  RECURVfiTURE 


Figure   1 2. 


Heidke   Skill    Scores    with    the    Value    of   the 
Predictand   from  ths    PSE   Equations    that   is 
Defined   tc  be    Recurvature   Varied    from   0.05    to 
0.95. 


smallest      variances      would   have      a      higher      HS3    due      to      the 
greater   probability      of   being   near   the     climatological    value 
(low    signal   regime)       than    m      the   tails    (high    signal   regime) 
of    the  predictand    distribution    (to   be   discussed   later). 

In  Table  XIII,  the  combined  P5E  method  produced  the 
highest  scores  relative  to  climatology,  except  at  the  36  h 
period  with  the  dependent  sez ,  where  the  PERS  method  was 
best.  An  indica-ion  of  the  clustering  of  predictands  about 
the      mean      appears      in      the      5^      h      independent,      set.  The 
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predictands  based  only  en  EOFS  have  the  highest  scores 
relative  to  climatology  at  this  time,  while  at  both  36h  and 
72  h  the  EOFS  were  inferior  to  the  other  two  methods.  P&E 
generally  showed  the  highest  skill  for  the  independent  set, 
with   the    exception   being      at   54    h  as   discussed   above. 

The  distribution  and  cumulative  frequency  of  the 
filtered  predictands  for  the  three  tests  are  given  in  Table 
XIV,  XV,  and  XVI.  An  unbiased  distribution  would  have  ^0% 
of  the  observations  in  each  ca-egory-  This  ideal  limit  is 
nearly  attained  at  the  36  h  time  period  for  the  dependent 
data  set.  However,  biases  become  more  evident  for  later 
time  periods.  The  t&E  test  generally  is  less  biased  than 
the  PERS  test.  Both  ?SS  and  PERS  are  Irss  biased  than  the 
EOFS  test,  which  again  is  an  indication  of  zhe  clusxering  of 
the    EOF    predictands   near    the    mean   value. 

TABLE   X  IV 

Fredicrand   Distribution   and    Cumulative    Percentage 
Distribution    for   36    h.      Cumulative   Percentage    is    in   Parens. 

Dependent  Set  Independent   Set 

36     h  PERS  EOFS  PSE  PERS  EOFS  ?SZ 

.0-  .1  36  (10)  38  (9)  41  (12) 

.  1-  .2  37  (21)  51  (21)  31    21 

.2-  .3  34(31)  70  (37)  42    33 

.3-  .4  35  (41)  69  (53)  40(44) 

.4-  .5  50(55)  47(66)  42(56) 

.5-  .6  32  (64)  43  (76)  31    65 

-6-  .7  33  (74)  33  (83  23  (71 

.7-  .8  32(83)  26(89)  22(73) 

.8-  .9  16  (87[  16  (93i  23(84) 

,9-1.0  44(100)  31(100)  55(100) 
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TABLE    XV 


Predictand   Distribution   and    Cumulative   Perce^Tage 
Distribution    for   54   h.      Cumulative   Percentage   is   m   Parens 


54    h 


Dependent  Set 
PERS  EOFS  PSE 


Independent   Set 
PERS  EOFS  P8E 


0- 

.  1 

1- 

.2 

2- 

.3 

3- 

.4 

4- 

.5 

5- 

.b 

6- 

,7 

7- 

.8 

8' 

.9 

9-1.0 


11 

20 
30 
45 
43 
39 
35 
26 
13 
9 


3 
12 
23) 

39 

55 

69) 
*8  2; 

92i 
.97' 
(100) 


111 

;i' 

! 

37 

2Q 

3 

42 

U3 

5 

42  ( 

.59 

14 

35 

72 

6 

27  ( 

82 

6 

14 

87 

1 

81 

90 

2 

28  ( 

[ID 

3) 

0 

TABI 

>  ill 

KVI 

14' 

21 
'33 
'65 
•79 

93) 
*95 

130) 

130 


6 

8 
10 
6 
6 
3 
0 
1 


^1 

'21 
[40 
*63 
[771 

91 

98] 
[98' 

100) 


]2 


9 
19) 

35i 
4  9 
70' 
81 
95] 
10  0) 
(100) 


predictand   Distribution   and    Cumulative    Pe^centaqe 
Distribution   for   72   h  .     Cumuli  ~  ' 


.ative   Percentage    is    in   Parens. 


72    h 

.0- 
.  1- 

.2- 

.  3- 
.4- 
.5- 
.6- 
,7- 
.8- 
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Dependent  Set 
PERS         EOFS  P&E 


Independent   Set 
PERS  EOFS  PSE 


.1 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
.0 
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The  Brier  P-3cor€s  for  the  three  individual  tests  and 
also  the  BPS  for  the  sample  climatology  and  for  chance  are 
given  in  Table  X7I I  which  indicates  that  each  of  the 
three  methods  is  superior  to  climatolDgy  and  c'^i-ixv^c^  at 
all  forecast  times.  This  is  encouraging  and  is  also  an 
indication  that  the  SCF  approach  of  describing  the  synoptic 
forcing        is        a        valid      approach        that         merits        further 
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TABLE  XVII 

F-Score  Statistics  for  36,  54  and  72  h  Dependent  and 
Independent  Data  Sets  After  Filtering  of  the  Predictand 
(A  Score  of  0.0  is  Perfecx) 


Depen 

dent  Set 

PEES 

EOFS 

P5E 

36  h   .3  08 

.348 

297 

climo 
chance 

.502 
.500 

54  h   .367 

.409 

340 

climo 
chance 

.500 

.500 

72  h   .389 

.451 

378 

climo 

chance 

.500 
.500 

Independent 

Set 

PERS 

EOFS 

P&E 

.390 

.454 

.490 
.500 

.383 

.401 

.393 

.504 
.500 

,368 

.443 

.441 

.500 
.500 

.417 

investigation.  Neither  PERS  Dr  EOFS  alone  shew  as  much 
skill  as  PSE  at  all  forecast  times.  This  is  true  for  both 
the  dependent  and  independent  data  sets.  The  difference 
between  the  two  data  samples  is  shown  by  the  large 
degradation  cf  the  BPS  at  36  h  for  rhe  independent  sample  in 
all  three  tests.  This  may  again  indicate  -hat  -he 
independent  sample  is  too  small.  The  scores  in  Table  XVII 
indicate  the  decrease  in  skill  as  time  increases.  I-  is 
believed  rhat  this  decrease  in  skill  with  respect  to  -ime 
could  be  lessened  to  some  degree  by  the  inclusion  of 
prognostic  EOF  fields.  This  was  indicated  by  Shaffer  (1982) 
as  one  of  the  largest  drawbacks  of  the  EOF  tropical  motion 
forecasting  procedure. 
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The  third  and  simplest  method  of  skill  analysis  is  the 
percent  correctly  forecast  (PCr)  method.  Table  XVIII  lists 
the  FCFs  obtained  for  each  test  for  both  the  dependent  and 
independent  data  sets.  This  table  again  illustrates  the 
superiority  of  the  P  5E  test  after  ths  initial  36  h  time 
period.  Notice  again  for  the  small  independent  data  se-c  at 
54  h,  the  EOF  method  is  superior  to  both  of  the  other 
methods  which  is  believed  to  be  due  to  the  clustering 
effect. 

TABLE  XVIII 

Percent    COiprfctly    Forecast    Eveats    from   th=    Regression    Method 
and  Defining   Recurvature  as   a    Prsdictand    Exceeding    the 
Climatology   Probabilities  of    Recurvai-ure 
(A    Score   of    100    is    Perfect! 

Dependent    set  Independent    Set 

PEES         EOFS  P€-E  PERS  EOFS  P&E 

36    h  77,7         76.0  76.3  72.9  65.3  72.9 

54    h  75.3         71.3  78.2  72.1  "76.7  72.1 

72    h  71,6         63.8  72.7  71.1  68.4  76.3 

This    simple   method     of  testing   skill    was      included    as    an 

attemp-c    to    compare      the   three    tested    procedures      against   the 

models   currently   being  used  by   J TWC   Guam    for   guidance.         The 

discussion   and      comments    pertaining      to   the      results    of      the 

model   comparisons    are   included   in    Appendix    A. 
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VI.   CONCLUSIONS  AND  JUSTIFICATION  FOR  CONTINUED 

RESEARCH 

Elsberry  (1979)  discusses  ths  inhsrent  problems  and 
advantages   of      dynamic  tropical    cyclone   models.  He    states 

that  the  lack  of  vertical  and  horizontal  resolution  in  the 
models,  as  well  as  a  poor  understanding  of  the  environmental 
processes  that  occur  in  the  vicinity  of  a  tropical  storm, 
are  limiting  the  progress  toward  more  accurate  forecasts. 
Although  nested  grid  models  (Harrison,  1981;  Harrison  and 
Fiorino,  1932)  have  been  used  to  improve  the  resolution, 
there  continue  to  be  problems  of  data  acquisition  on  the 
scales  of  interest.  This  problem  has  no  solution  in  the 
immediate  or  extended  time  frames,  which  necessitates  the 
need    for    other   reliable    and  less    costly    methods. 

The  application  of  a  linear  regression  scheme  t.o  a 
statistical  model  using  a  combination  of  persistence-related 
parameters  and  an  EOF  representation  of  the  synoptic  forcing 
was  completed.  The  combination  of  these  two  types  of 
parameters  have  shown  to  be  superior  for  forecasting 
tropical    cyclone  recurvature   as    defined. 
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The  main  attributes  of  th=  EOF  approach  to  foracasting 
tropical  cyclone  recurvature  and  motion  (Shaffer, 1982)  must 
again   be   stressed: 

•  The    method    is   extremely   inexpensive; 

•  The  forecasts  can  be  done  on  a  haadheld  programmable 
calculator; 

•  The  method  has  been  shown  to  be  more  accurate  for 
forecasxing  recurvature  than  persistence  and  synoptic 
parameters    alone;  and 

•  The  method  has  shown  to  give  superior  forecasts  when 
compared  to  the  official  JTWC  "forecasts  using  the 
analysis   of   the    500    mb  fields    (Shaffsr,     1982). 

It  is  believed  that  the  inclusion  of  prognostic  synoptic 
fields,  as  suggested  by  Shaffer  (1982),  would  improve  the 
percentage  of  the  explained  variance  of  the  predictand. 
This  would  also  be  reflected  as  giving  batter  results  for 
the    recurvature   study   described    here. 

An  independent  operational  evaluation  of  the  schemes 
presented  here  needs  to  be  done  to  determine  the  validity  of 
the    method. 

The  definition  of  recurvature  needs  to  be  modified,  k 
tropical        storm        moving        314°  is        considered        as        a 

non-recurvatur e  case  where  a  storm  moving  3 1 6  <^  is  coasidsred 
a  recurvature  case.  Probabilities  of  recurvature  based  on 
the    departure    of   the    tropical      cyclone    crDin    a    forecast    track 
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or  a  straight  line  could  be  utilized  instead  of  using  a 
binary  definition  as    done    in   this  rssearch. 

It  is  understood  that  the  application  of  new  and 
superior  methods  of  forecasting  tropical  cyclone  recurvanure 
and    notion   should   be    questioned.  It   is   believed   rhat    this 

change  of  methodology  (ie,  using  an  EOF  representation  of 
the  synoptic  fields),  if  adopted,  will  isprove  the  image  of 
the  military  weather  services  by  being  a  positive  pioneering 
effort  in    an   area   that  has   been    stagnant    for    many   years. 
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APPENDIX    A 
MODEL   COMPARISONS 


A.  PURPOSE 

The   purpose     of   this      section   is      to    supply      the    Western 

Pacific  tropical   cyclone    forecaster      with    information    on   the 

various  guidance   techniques   for    forecasting    tropical    cyclone 

motion,    and,    in    particular,   their   reliability   in    forecasting 

recurvature.  The      definition      of      recurvailure      (Leftwich, 

1978,1982)       is    similar      as  described    earlier    except      for    the 

northeastward        motion     during      the   24   h    (vice    12)      prior    to 

the    end    of    t.he    forecast    period.         This   modification   was   done 

to    accomoda-e  the' forecast      warning    verification    times    which 

are    24,    48    and   72    h. 

"Recurvature  is  a  nex  displacement  northward  of  315° 
during  the  forecast  period  or  the  attaining  of 
northeastward  motion  during  the  24  hours  prior  to  tae  end 
of   ths    forecast    period. 

B.  THE    DATA 

The  data  analyzed  were  assembled  by  Dr.  Ted  Tsui  of  the 
Naval  Environmental  Prediction  Research  Facility  (NEPR?) 
during  the  1979-1981  Northwest  Pacific  typhoon  seasons.  The 
data      were   the      best      track      information   from      the      trooical 
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disturbance  stage  through  the  entire  lifs  cycle  of  each  of 
the  28  storms  of  1979  and  1980  and  the  29  storms  of  1981. 
Also  included  were  the  24,  48  and  72  h  positions  of  12 
tropical  cyclone  forecast  aids  available  at  JTWC.  Seven  of 
these  aids  and  the  Joint  Typhoon  Warning  Center  official 
forecasts  are  evaluated  below  with  respect  to  the 
recurvature   decision. 

C.        MODEL    DESCRIPTIONS   COMPARISON    TESTS 

The  following  paragraphs  are  das^riptions  of  the 
objective  forecast  aides  currently  being  used  at  JTWC.  The 
descriptions  have  been  taken  from  the  199  1  Annual  Tropical 
Cyclone  Report.  All  JTWC  objective  techniques,  except  -chose 
involving  dynamic  models,  are  executed  using  operational 
data  available  prior  to  warning  times  (OOGMT  and  12g:iT). 
These  dynamically-based  objective  aids  are  usually  received 
at    JTWC    within   2   to   5    hours   of  a    specific    warning   time. 

An  analogue  method,  abbr=via-3d  TOPS,  provides  three 
movement  and  intensity  forecasts.  These  three  forecasts  are 
for  analogues  of  st raighx-movi ng  cyclones,  of  rscurving 
cyclones  and  a  combination  of  analogues  including  s-raight, 
recurving,  and  all  other  cyclones  that  do  not  meet,  the 
specific   cirteria   of    those  two  categories. 
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An  extrapolation  technique,  abbreviated  XTRA,  provides  a 
track  connecting  the  12-h  old  preliminary  best  track 
position  and  the  curi;ent  warning  position,  which  is  then 
extrapolated  to  24  and  43  h.  HPAC,  which  is  a  combination 
of  XTRA  and  a  cl im at o logical  aid,  provides  a  forecast  of  the 
track  and  intensity  changes  based  on  the  initial  position  of 
the  system.  The  c lima tologiz al  data  base  includes  the 
period  from  1945  to  1973.  The  ( BPAC)  forecasts  are  based  on 
blending  the  past  motion  of  the  rropical  cyclone  with  the 
climatology  forecast  positions.  This  program  requires  only 
a  Texas  Instrument  (TI-59)  calculator  to  generate  12  to  72  h 
forecast  positions.  The  blending  routine  gives  less  weight 
to    persistence   at   each  succeeding   forecast    interval. 

An  updated  version  of  the  HATRACK/tiOHATT  steering 
program  (now  called  GfCLOPS  and  designated  as  CY50  for  the 
500  mb  level)  can  provide  steering  forecasts  at  niany 
different  levels.  The  program  can  be  run  in  the  modified 
(includes  a  12-h  persistence  bias)  or  anmodified  versions 
applied      to      either      analysis    or      prognostic      fields.  The 

program  advects  a  point  vortex  on  a  pre-selected  analysis 
and/cr  smoothed  prognostic  field  at  designated  levels  in  o-h 
time   steps      through  72      h.        In      the    modified      version,      the 


64 


program  us€s  the  previous  12-h  history  position  -o  compute 
the  12-h  forecast  error  and  applies  a  bias  correction  to  the 
forecast  positions.  In  1981,  only  the  modified  version,  in 
the    prognostic   mode   for   the  500    mb   level    was    verified. 

.The  dynamic  Tropical  Cyclone  Model  (TCMO)  is  a  coarse 
mesh  (220  km)  primitive  equation  model.  The  digitized 
cyclone  warning  position  is  bogused  into  the  850  mb  wind  and 
temperature  fields  of  the  F NOC  global  band  analysis. 
Hemispheric  model  forecast  data  are  used  on  the  boundaries. 
Two  versions  are  currently  run:  One  runs  from  forecast 
fields  and  is  available  on  an  automatic  request  basis,  and 
the  other  is  initialized  on  the  analysis  fields  not^d 
earlier.  The  model  based  on  the  analysis  fields  was  tested 
in    this   study. 

The   NTCM   or   "nested"    primitive   equation    tropical   cyclone 
model   is      initialized   on      FNOC    1  2-h      forecast,    fields.  The 

model  covers  a  limited,  but  relocatable,  tropical  domain 
with  three  layers  in  the  vertical.  The  finer  scale,  or 
"nested"  grid,  covers  a  1200  k.  m  square  area  with  a  41  km 
grid  spacing  and  moves  to  keep  a  850  mb  vortex  in  its 
center.  This  model  is  integrated  with  a  coarse  channel 
model   grid    with   a   grid  spacing   of    20  5    km    over   a    6400    Dy    4700 
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km    domain.         Once   initialized,    the    model   runs   independent:   of 
the    remaining   FNOC   forecast   fields.  The   NTCM    is   available 

via      automatic   request      for  00      and    12      GMT    forecast      fields 
only. 

D.        SCREENING    PEOCEDUHE    FOE   MODEL    TESTING 

Several  diffarert  requirements  were  taken  into 
consideration  for  inclusion  in  this  comparison  study  that 
were  not  considered  for  the  EOF  analysis  scheme.  Inclusion 
was   tased   on: 

•  Using      the      'best      track'       past      zonal      and      meridional 
displacements    (after      the   initial      12    h      of   the      track) 
vice    the      past    warning      displacements    used      in   the      liOF 
study,        to        test      for        eastward      displacement         (not 
consiaered   case)  ; 

•  Pseudo-ir.  depende  nee    is  again   assumed      by   requiring    time 
separations' of    3  6  h    within    the   same    storm; 

•  A  30    h    past    history    was    not    required; 

•  All   tropical   cyclones   wer^      included,      rather  than    just 
tropicaj.   storms    and    typhoons;    and 

•  All    storms    north  of    the    equator,      vi^s    10°    north,      were 
included. 

Table      XIX    lists      the      24,       48  and      72    h      recurvature 

probabilities    for   the     three    year   sample.       When      compared    to 

the    corresponding    values    of   the    sample   climatology  listed    in 

Table     II,      the      difference  between      tha    inclusion      criteria 

become   evident.         The    inclusion      of    the   tropical    depressions 
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in  the  area  south  of  10°  north  is  the  largest  contributor  to 
the  lower  values  of  rscurvaturs  probability,  A  high  annual 
variability   is      indicated   in   this     small   sample.  There    is 

generally  a  slow  increase  of  probability  Df  recurvature  with 
time   as   was   the   case    in   Table   II . 

TABLE   XIX 

Nuffibep  9f   Recurving    (rev),    Non-recurvina    /nrcv),    and 
Probability  of   Recurvature    (prob)     for    Tropical  Cyclone 
Seasons    1979,    1980    and    1981. 

1979  1980  1981 

tau      rev      nrcv  prob  rev  nrcv  prob  rev  nrcv  prob 

24         39           94  .293  36  78  .316  26  73  .263 

48         41            84  .328  35  55  .350  23  63  .267 

72         25           81  .236  29  53  .354  21  52  .268 

The  HSS  for  the  three  separate  years  are  listed  in  Table 
XX.  The  HSS  in  this  data  set  remains  as  described  earlier 
with  a  slight  modification.  This  difference  is  tnat  the 
verification  of  the  model's  pes  dieted  recurvature  is  based 
entirely  upon  the  predict.ands  matching  the  actual 
recurvature  criteria  as  defined.  There  is  no  climatological 
threshhold  that  is  set  for  in::lusion  as  a  recurvature  or 
non-recurvature    as   was  the   case    for  the    earlier   data    sample. 

Again,  the  most  noticable  feature  of  the  table  is  the 
general      decrease      of      skill    with      respect      to      time.  The 
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exception  to  this  generality  is  the  3PAC  model  performance 
for  1979  and  the  TCMO  model  for  1981  where  the  HSS  increased 
with  time.  An   interesting      comparison   i^an      be   shown      when 

comparing  JTWC  with  the  TCMO  forecast  aide  for  1981.  There 
is  a  definite  correlation  between  TCMO  ani  JTWC.  This  is  an 
indication  of  JTWC  fcllowing  the  guidance  of  this  dynamic 
model   quire   closely   in  the   extended  tiae    frames. 

Table  XXI  lists  the  BPS  for  the  three  time  periods  of 
interest.  The  BPS  forecast  predictands  in  this  case  are  not 
on  the  interval  between  1.0  and  0.0,  but  are  identically 
equal  to  1.0  or  0.0.  Any  knowledge  of  a  comparison  of  BPS 
results  using  a  continuous,  predictand  as  compared  to  a 
binary  predictand  scheme  is  available.  Comparison  of  the 
common  time  period  of  72  h  between  the  -wo  en-irely 
independent  data  sets  of  diff=rem:  sizes  there  fore  has 
little  meaning.  The  BPS  of  tha  TCMO  model  in  1981  and  BPAC 
for  1979  are  the  only  models  tha*  approach  the  BPS  for  the 
combination  of  EOFS  and  Persistence  in  the  earlier  study.  A 
general  trend  toward  larger  BPS  are  prevalent  for  all  years 
for    most    of   the    forecast    aids. 

Comparison  of  the  PCF  values  indicated  in  Table  XXII, 
and    the    HSS   reveals   the   similarity      of   -he    two    skill    scores. 
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TABL2     XX 

Heidke    Skill    Sco^e    for  Models    for    1979.     1980,    and    1931 

(A  Score  Of    1.9  is   Perfect,    and  a    (*)  indicares   too 

Few   Cases  to   be  Ranlced) 

Model  results   for  1979 

HPAC         TCMO  JTHC  TOPS  CY50  B?AC  XTRA  NTCM 

2a    h      .1*49         .521  -424  .296  .388  .403  .304  .167 

cases      421              95  427  247           3  04  345  425  12 

rank        (2)            (1)  (3)  (7)            (5)  (4)            (6)  (*) 

48    h      .457         .506  .355  .230  .426  .568  .235  .667 

cases      368              81  363  228           265  296           371  12 

rank         (3)            (2)  (5)  (7)            (4)  (1)            (6)  (*) 

72    h       .381         .391  .347  .266  .095  .586  .072  .657 

cases      252              69  288  199           210  251           304  12 

rank         (3)            (2)  (4)  (5)            (6)  (1)             (7)  (*) 

Model  results    for  19  80 

HPAC         TCMO  JTWC  T0P5  CYSO  BPAC  XTEA  NTCM 

24    h       .518         .621  .424  .296  .388  .403  .304  .167 

cases      419              73  412  353           364  361           373  62 

rank        (2)            (1)  (3)  (7)            (5)  (4)            (6)  (8) 

48    h      .313         .406  .418  .461  .411  .453  .301  .387 

cases      352             64  347  3  04           309  300           306  62 

rank         (7)            (5)  (3)  (1)            (4)  (2)             (8)  (6) 

72    h       .004         .378  .  C89  ,112  .273  .217  .164  .309 

cases      261              45  257  245           245  240           225  55 

rank         (8)            (1)  (7)  (6)            (3)  (4)            (5)  (2) 

Model  results   for  1981 

HPAC         TCMO  JTWC  T0P5  CY50  BPAC  XTRA  MTCM 

24    h       .462         .457  .542  .466  .567  ,480  .483  .333 

cases      383            129  384  348           282  377           383  156 

rank        (6)            (7)  (2)  (5)            (1)  (4)            (3)  (8) 

48    h      .298         .574  .515  .333  .439  .388  .255  .289 

cases      325            108  309  297           235  320           325  135 

rank         (6)            (1)  (2)  (5)             (3)  (4)            (8)  (7) 

72    h      .013         .647  .540  .021  .304  .063  ,149  .196 

cases      235              68  268  235            181  254           235  107 

rank         (8)            (1)  (2)  (7)            (3)  (6)             (5)  (4) 

The    rankings   of   these   two  methois  are      th=    same  for    51  of    69 

comparisons. 
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TABLE   XXI 


Brie 
1 


579/1 

a 


core   Staxistics  for    2U,    U8    api    72    h    Models    for 
80   and    1981     (A  Scora    of   0.0    is    Perfect,    and 
(*)    Indicates   too  Few    Cases   to   be    Ranked) 


Model  results   for    1979 


24    h 

cases 

rank 

US    h 
cases 

rank 

72    h 

cases 

rank 


HP  AC 

.570 
421 
(2) 

.543 
368 
(3) 

.6  19 
252 
(3) 


TCMO 

.379 

95 

(1) 

.494 

81 

(2) 

.609 
69 
(2) 


HP  AC 
24  h  ,482 
cases  4 19 
rank        (5) 


48    h 
cases 

rank 


688 
352 
(7) 


72  h  1.004 
cases  261 
rank        (8) 


TCMO 
.493 
73 
(6) 

.  594 

64 

(5) 

.622 
45 
(1) 


HP  AC 
24  h  ,538 
cases  383 
rank         (6) 


48    h 
cases 

rank 

72    h 
cases 

rank 


.702 
325 
(6) 

,9  87 
235 
(8) 


:cMo 
543 

129 
(7) 

426 
108 

(1) 

,35  3 

68 

(1) 


JTWC 

.576 
427 
(3) 

.  645 
363 
(5) 

.  653 
288 


Model  results   for    19  80 


JIWC 
.  461 
412 
(3) 

.  582 
347 

(3) 

.  911 

257 
(7) 


Model  results   for   19  31 


JIWC 
.  458 
384 
(2) 

.  485 
309 
(2) 

.  460 
268 
(2) 


T0P5 

.704 
247 
(7) 

.702 
223 
(6) 

.734 
199 
(5) 


T0P5 

.465 

353 

(^) 

.539 
3  04 

(1) 

.873 
245 
(6) 


T0P5 

.534 

348 

(5) 

.667 
297 
(5) 

.9  79 
235 
(7) 


CY50 

.612 
3  04 
(5) 

.574 
265 
(4) 

.905 
210 
(6) 


CY50 

.440 

364 

(1) 

.589 
309 

(^) 

,727 
245 
(3) 


CY50 

,433 

282 

(1) 

.511 
235 
(3) 

.696 
181 
(3) 


3PAC 

,597 
345 
C*) 

.432 
296 

(1) 

.414 
251 
(1) 


BPAC 
.454 
361 
(2) 

.547 
300 
(2) 

,783 
240 
(4) 


3PAC 

,  520 

377 

(4) 

.612 

320 

C^) 

,937 
254 
(6) 


XTEA 

.6  96 
425 
(6) 

.765 
371 

(7) 

.928 
304 
(7) 


XTRA 
.536 
373 
(7) 

.699 
306 

.836 
225 


XT  RA 

.517 

383 

(3) 

,745 
325 
(8) 

,902 
235 
(5) 


NTCM 

.333 
12 

(*) 

,333 
12 

(*) 

,333 

12 

(*) 


NTCM 
.581 
62 
(8) 

.613 
62 
(6) 

.691 

55 

(2) 


NTCM 
.667 
156 
(8) 

,711 
135 
(7) 

.804 
107 
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TABLE    XXII 


Pe 


rcent   Correct   Forecast   Statistics   for    Models   for    1979, 
1980    and    1981       (A   Score   Of    100   is    Perfectj    and  a    (*) 
Indicates   t  co    Few  Cases   to    be    Rankea) 


Model  results    for    1979 


HP  AC 


TCMO 


JTWC         TOPS 


CY50 


BPAC 


XTRA 


NTCM 


2U    h 
cases 

rank 

71.5 

421 
(2) 

81.  1 

95 

(1) 

7  1.2 

427 

(3) 

2  47 
(7) 

69.4 
3  04 
(5) 

70.  1 

345 

65.2 
425 
(6) 

58.  3 

12 

(*) 

48    h 

cases 

rank 

72.8 
368 
(3) 

75.  3 
81 
(2) 

6  7.8 
363 
(5) 

64.9 
2  28 
(6) 

71  .3 

265 

78.4 
296 
(1) 

61.7 
371 
(7) 

83.  3 

12 

(*) 

72    h 
cases 

rank 

(3) 

69.6 

69 

(2) 

6  7.4 
288 

63.3 
199 
(5) 

54.8 
210 
(6) 

79.3 
251 
(1) 

53.6 
304 
(7) 

83.3 
12 

(*) 

Model  results   for 

19  SO 

24    h 

cases 
rank 

HP  AC 
75.9 
419 
(5) 

TCMO 

75.3 

73 

(6) 

JTWC 
7  6.9 
412 
(3) 

TOPS 

76.8 
353 
(4) 

CY50 
78.0 
364 
(1) 

BPAC 
77.3 
361 
(2) 

XTRA 

73.2 

373 

(7) 

NTCM 

71.0 

62 

(8) 

43    h 

cases 

rank 

65.6 
352 
(7) 

70.3 
64 

7  0.9 
347 
(2) 

64.9 
3  04 
(8) 

70.6 
309 
(3) 

72.7 
300 
(1) 

69.3 
306 
(6) 

69.4 

62 

(5) 

72    h 

cases 

rank 

49.8 
261 
(8) 

68.9 
45 

(1) 

5  4.5 
257 
(7) 

63.3 
245 
(5) 

63.7 
245 

6  0.8 
240 
(6) 

64.0 
225 
(3) 

65.4 

55 

(2) 

Model  resul- 

ts   for 

19  81 

- 

24    h 
cases 

rank 

HP  AC 
73.  1 
383 
(6) 

TCMO 
72.9 
129 
(7) 

JTWC 
7  7.1 
384 
(2) 

TOPS 

73.3 

3  43 

(5) 

C£50 
78.4 
282 
(1) 

BPAC 

7  4.0 

377 

XTRA 
74.2 
383 
(3) 

NTCM 
66 ,  7 
156 
(8) 

48    h 
cases 

rank 

64.9 
325 
(6) 

78.7 
103 
(1) 

7  5.7 
309 
(2) 

66.7 

297 
(5) 

74.5 
235 
(3) 

69.4 
320 

62.8 
3  25 
(8) 

64.  4 
135 

72    h 
cases 

rank 

50.6 
235 
(3) 

83.4 

68 

(1) 

7  7.0 
268 
(2) 

51.1 
235 
(7) 

65.2 
181 
(3) 

53.1 
254 
(6) 

54.9 
235 
(5) 

59.8 
107 

E. 

SUMMARY 

Seven        forecast         aids      used        by        Joim:  Typhocri 

Warning   Canzer       (JTWC)      Guam,         were    tsstad      for   ability      to 
forecast   tropical   cyclone   recurvature   at   24,      48        and   72    h. 
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Using  the  skill  scores  described  in  Chapter  V  for  the 
1979-81  typhoon  season,  model  performance  indicates  a  large 
annual  variation   of   skill   between   the   modsls. 
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